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PREFACE 


The  advent  of  slender  high-speed  re-entry  bodies  and  particularly 
slender  lift-producing  surfaces  re-emphasizes  the  importance  of  phe¬ 
nomena  occurring  near  the  leading  edge.  Previous  work  at  The  RAND 
Corporation  already  considered  such  phenomena  as  the  coupling  of  hy¬ 
personic  viscous  induced  pressure  and  mass  addition  at  the  surface  — 
typically  a  leading-edge  problem  showing  measurable  effects  on  the 
aerodynamic  force  field.  However,  that  analysis  was  based  on  an  Ide¬ 
alized  model  of  the  flow,  namely,  on  classical  similarity  solutions 
for  binary  boundary  layers. 

This  Memorandum  explores  the  problem  of  the  coupling  between  the 
phase  change  of  a  sublimating  surface  and  the  flow  near  the  leading 
edge.  The  primary  value  of  the  results  is  the  establishment  of  the 
nature  of  this  coupling  and  of  ^ts  importance.  The  results  contribute, 
for  example,  to  the  decision  regarding  the  conditions  under  which  it 
i s  realistic  to  proceed  with  studies  of  hypersonic  induced  pressure 
interaction  with  surface  blowing  without  simultaneously  including  the 
coupling  of  the  blowing  rate  with  the  flow  field. 
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SUMMARY 


The  behavior  of  the  laminar  binary  boundary  layer  with  blowing 
is  coupled  with  the  kinetic-evaporation-rate  law  of  the  surface  mate¬ 
rial.  It  is  found  that  the  solution  of  this  combined  problem  exhibits 
an  asymptotic  behavior  at  large  Reynolds  number  which  is  identified 
with  the  usual  near-equilibrium  solution  in  which  the  vapor  net'r  the 
wall  is  verv  nearly  in  thermodynamic  phase-change  equilibrium.  How¬ 
ever,  the  near-equilibrium  solution  is  invalid  in  a  region  near  the 
origin  of  the  boundary  layer,  which  is  characterized  by  a  length  formed 
with  physical  parameters  describing  the  basic  flow  and  the  surface 
properties.  In  this  region,  which  is  treated  here  approximately,  the 
blowing  parameter  decreases  to  zero,  the  wall  temperature  increases, 
and  the  sublimation  rate  tends  to  an  upper  limit  as  the  origin  of  the 
boundary  layer  is  approached. 

The  principal  result  of  the  analysis  is  the  estimation  of  the 
length  of  the  region  of  transition  to  the  near-equilibrium  solution. 
This  length  is  an  independent  characteristic  of  the  problem  and  does 
not  scale  as  do  the  boundary-layer  properties  (that  is,  with  the  Rey¬ 
nolds  number).  This  implies  that  while  the  transitional  length  on 
typical  re-entry  vehicles  is  small  and  probably  negligible  unless  the 
leading  edge  is  very  sharp,  practical  wind-tunnel  test  models  can 
easily  be  affecieii  over  most  of  their  chord. 
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I .  INTRODUCTION 


The  complete  solution  of  a  flow  field  over  a  subliming  surface 
represents  an  equilibrium  among  the  rate  of  heat  tnobfer  to  the  wail, 
the  rate  of  phase  change  of  the  surface  material,  and  the  rate  of  dif¬ 
fusion  of  the  vapor  evolved  at  the  wall  through  the  boundary  layer. 
These  conditions  determine  the  surface  temperature  and  the  blowing 
rate. 

The  problem  is  usually  treated  as  follows:  The  binary  boundary- 
layer  equations  are  solved  with  the  wall  temperature  and  the  blowing 
rate  treated  as  independent  boundary  conditions.  The  solution  yields 
the  temperature  gradient,  that  is,  the  heat-transfer  rate,  as  a  para¬ 
metric  function  of  the  wall  temperature  and  the  blowing  rate.  Since 
the  blowing  rate  and  the  heat  transfer  are  related  by  the  effective 
heat  of  sublimation  of  the  material,  one  more  relation  is  needed  to 
fix  uniquely  the  two  free  parameters  (wall  temperature  and  blowing 
rate).  At  this  p^int  we  make  the  approximation  that  the  phase  change 
occurs  at  thermodynamic  equilibrium;  that  is,  the  wall  temperature  is 
the  phase-equilibrium  temperature  at  the  existing  partial  pressure  of 
the  vapor  near  the  wall  (which  is  known  from  the  solution  of  the  bina¬ 
ry  boundary- layer  equations).  This  relation  suffices  to  complete  the 
formulation. 

The  assumption  of  thermodynamic  phase  equilibrium  at  the  surface 
is  conceptually  incorrect  because  at  equilibrium  the  net  mass  transfer 
between  phases  Is  zero.  Therefore,  a  more  complete  kinetic  relation 
among  the  tempiraturt .  concentration,  and  rate  of  phase  change  is  ac¬ 
tually  needed.  This  has  been  discussed  in  a  number  of  papers, 
hu:  no  solution  including  such  s  hio*tic  surface- evaporation  condition 
has  been  obtained,  nor  have  the  implication  of  this  phenomenon  been 
fully  explored. 

The  nesr-equilibrlum  solution  is  valid  at  sufficiently  high  Key- 
no  Id*  numbers,  which  cm  be  illustrated  by  saying  that  when  rhe 
Reynolds  number  is  high,  the  lapedance  of  the  boundary  layer  to  dif¬ 
fusion  of  vapor  is  very  much  higher  than  the  lapedsnce  to  surface 
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(4) 

phase  change,  which  is  in  series  with  it.  The  latter  is  then  neg¬ 
ligible,  and  the  coupling  between  the  flow  and  the  surface-evaporation 
phenomenon  is  diffusion  limited.  This  statement  leads  to  the  question, 
What  exactly  are  the  lower  limits  of  validity  of  this  approximation? 

If  the  surface-evaporation-rate  law  is  considered,  the  results 
of  the  near-equilibrium  approximation  must  be  interpreted  as  saying 
that  the  difference  between  the  actual  wall  temperatuu  and  the  thermo¬ 
dynamic-equilibrium  temperature  of  the  vapor  is  small  everywhere.  The 
actual  wall  temperature  is  determined  by  the  rate  equation  so  as  to 
supply  the  mass  flux  from  the  solid  to  the  vapor  phase.  The  equilib¬ 
rium  temperature  is  determined  by  the  partial  pressure  whicn  must  ex¬ 
ist  near  the  wall  to  drive  the  mass  flux  across  the  boundary  layer  by 
pressure  diffusion.  Smallness  is  measured  in  comparison  with  the 
temperature  difference  driving  the  heat  transfer;  that  is,  .he  differ¬ 
ence  between  adiabatic  recovery  and  wall  temperature. 

The  well-known  result  of  near-equilibrium  solutions  for  self¬ 
similar  flows,  for  instance,  is  that  the  surface  (equilibrium)  temper¬ 
ature  is  constant,  while  the  sublimation  rate  varies  as  the  inverse 
square  root  of  the  Reynolds  number.  We  may  consider  the  Reynolds  num¬ 
ber  as  a  unique  measure  of  the  distance  from  the  origin  of  the  boundary 
layer  if  the  free  stream  is  fixed  and  the  wail  temperature  is  constant 
Rut  if  the  mass  flux  increases  towards  the  origin  of  the  boundary  layer, 
then  the  wail  temperature  must  increase  according  to  the  sublimation- 
rate  law.  THo  vail  temperature  cannot  he  constant,  and  therefore  the 
Reynolds  number  cannot  be  a  unique  measure  of  distance.  Also,  the 
near-equilibrium  approximation  must  break  down  at  some  distance  to¬ 
wards  the  origin. 

The  properties  of  the  boundary  layer  scale  in  terms  of  two  param¬ 
eters,  the  Reynolds  number  and  the  wall  temperature  more  correctly, 
some  ratio  of  free-strea*  temperature  to  wall  temperature)  Having 
recognised  that  the  wall  temperature  is  coupled  to  the  problem  through 
an  Independent  rate  eqvition,  one  concludes  that  the  introduction  of 
the  kinetics  of  surfat-  phase  change  brings  in  *  new  scale  parameter 
This  scale  has  the  fe  a  of  a  distance  from  the  origin  of  the  boundary 
layer  for  fixed  free  stream  conditions  and  given  thermodynamic  phase- 
change  conataata. 
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The  condition  for  validity  of  the  near-equilibrium  approximation 
can  now  be  stated  as  follows:  The  approximation  is  valid  when  the 
Reynolds  number  is  much  higher  than  the  Reynolds  number  formed  with 
the  leading-edge  scale  described  above.  Moreover,  the  Illustrative 
argument  based  eo  the  relative  magnitude  of  the  impedance  to  mass 
transfer  is  seen  to  be  incomplete;  it  does  not  consider  the  nonlinear 
coupling  between  the  "impedance”  and  the  driving  potential,  for  in¬ 
stance,  it  implies  that  when  the  Reynolds  number  is  verv  low  the 
boundary- layer  transfer  impedance  is  low  and  the  coupling  between  flow 
and  sublimation  is  dominated  by  the  surface  impedance  (rate  limited) 
This  is  not  true  if  the  Reynolds  number  is  low  by  virtue  of  low  pres¬ 
sure  at  a  given  distance  from  the  origin.  It  is  only  true  if  the 
Reynolds  number  is  low  by  virtue  of  small  distances  from  the  leading 
edge . 

The  purpose  of  this  study  is  to  investigate  the  character  of  the 
rate-limited  sublimation  problem  and,  in  particular,  to  determine  the 
scale  of  this  region.  It  is  proper  to  make  two  remarks  at  this  point. 
First,  the  analysis  is  based  entirely  on  continuum  boundary-layer 
concepts  This  is  equivalent  to  saying  that  the  scale  of  the  rate- 
limited  region  must  be  larger  than  some  minimum  distance  required  for 
validity  of  continuum  boundary-layer  concepts  for  the  analysis  to  be 
meaningful.  This  situation  turns  out  to  be  practically  possible. 
Second,  in  regard  to  practical  flows  over  somewhat  blunted  bodies, 
the  statement  "origin  of  the  boundary  lay**r”  or»t  be  interpreted  a* 
meaning  a  virtual  origin  from  which  the  boundary  layer  would  start  to 
attain  a  thickness  and  profile  it  has  at  the  point  under  cons ' derat  ion 
This  implies  that  the  distance  to  this  virtual  origin  mist  be  larger 
than  the  radius  at  th"  blunted  nose. 

The  current  trend  towards  finer  re-entry  shape*  and  the  concern 

with  problems  such  as  the  interaction  between  ablative  mass  addition 

(b) 

and  hypersonic  viscous  induced  pressure  --typically  a  leading-edge 
phenomenon --tends  to  bring  t be  transitional  sublimation  regime  into 
the  realm  of  practical  problems 

Finally,  the  present  solutions  for  the  rate-limited  sublimation 
region  involve  a  number  of  approximations  in  the  treatment  of  the 
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boundary  layer  and  should  be  interpreted  mainly  as  a  study  of  whether 
or  net  the  problem  is  sufficiently  significant  to  deserve  a  more  rig¬ 
orous  and  much  more  difficult  analysis. 
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II.  FORMULATION 


The  problem  requires  the  simultaneous  solution  of  both  the  bound¬ 
ary  layer  and  surface  phase- change  rate  equations.  Phenomenologically, 
these  fall  into  lour  groups  describing,  respectively,  the  energy  and 
the  mass-transfer  properties  of  the  boundary  layer  and  the  surface. 


BOUNDARY-LAYER  HF.AT  TRANSFER 

The  first  is  a  solution  of  the  classical  binary  boundary-layer 
equations  which  we  consider  to  be  uncoupled  from  the  mass-diffusion 
equation, ^  implying  a  Lewis  number  approximately  equal  to  one.  We 
do  not  consider  chemical  reactions  between  the  interdif fusing  species. 
The  solution  has  the  form 


where  B  is  a  normalized  blowing  parameter 


B 


m 


p  u  C, 
Ke  e  h 


(1) 


(2) 


The  subscript  o  indicates  reference  conditions  (nonablating  surface), 

and  the  Stanton  number  C,  is  defined  on  the  basis  of  the  adiabatic 

n 

recovery  enthalpy  (temperature): 


'h  p  u  (h  -  h  ) 
e  e  aw  w 


p  u  c  (T  -  T  ) 
Ke  e  p  aw  w 


(3) 


Equation  (1)  formally  includes  the  influence  of  the  pressure  gradient, 
the  temperature  gradient,  the  variation  of  mean  molecular  weight,  and 
gas  properties  through  the  boundary  layer.  This  is  reflected  in  the 
expression  for  the  function  F.  Practically,  only  similarity  solutions 
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are  known,  and  for  these  the  well-known  linear  approximation  for  F  in 
terms  of  B  holds  to  values  of  B  on  the  order  of  0.3.  The  numeric^1 
value  of  the  proportionality  constant  and  its  dependence  on  Prandtl 
number,  mass  ratio,  and  pressure  gradient  (within  the  similar-flows 
family)  has  been  discussed  exhaustively  in  the  literature . ^  ^ 

We  propose  to  write  an  approximation  to  Eq.  (1): 


(4) 


Equation  (4)  is  co  be  considered  as  a  formal  approximate  expression  of 

the  behavior  of  the  Stanton  number  for  B  sufficiently  small,  with  u  an 

unspecified  function  of  the  indicated  parameters.  The  dependence  on 

(9) 

the  Prandtl  number  is  suggested  by  the  similar  solutions. 

It  will  be  seen  a  posteriori  that  E  goes  from  0  when  Re^  *  0  to 

an  asymptotic  value  when  Re^  -♦  In  the  initial  region  the  wall- 

temperature  gradient  is  very  strong  and  local  similarity  not  valid; 

but  Eq.  (4)  is  then  defencable  as  the  first  term  of  an  expansion  about 

* 

B  *  0.  In  the  region  of  asymptotic  approach,  gradients  are  small  and, 
provided  B^  <0.3,  Eq.  (4)  holds  as  a  "local  similarity"  solution.  It 
is  difficult  to  imagine  a  situation  in  which  essential  errors  in  the 
trends  exhibited  in  the  transition  region  would  result  from  the  use  of 
Eq.  (4). 

To  the  same  degree  of  approximation  the  recovery  factor  for  the 
binary  boundary  layer  is  expressed  in  terms  of  the  recovery  factor  for 
the  reference  flow: 


The  external  pressure  gradient  is  an  independent  parameter.  It 
is  assumed  to  be  such  that  use  of  the  "local-similarity"  concept  can 
be  justified  in  regard  to  it. 
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r 


(5) 


The  numerical  values  of  g  *  g(Pr,  dp/dx,  dT  /dx,  M./M  )  are  less  well 
defined  from  available  solutions  even  for  self-similar  flows.  It  will 
be  seen  later,  however,  that  the  proportionality  constants  at  and  B  ap¬ 
pear  in  the  solution  only  as  a  ratio.  Thus,  if  they  both  vary  in  the 
same  way  with  M  /Mi  Pr,  dp/dx,  and  dT  /dx  to  a  first  approximation, 
the  effect  of  this  dependence  is  minimized. 

THE  HEAT  BALANCE  AT  THE  WALL 

The  thermal  coupling  between  the  flow  field  and  the  boundary  is 
specified  by 


q 


m  X 
s  s 


(6) 


where  X^  is  the  heat  of  sublimation  of  the  wall  material.  This  approx¬ 
imation  neglects  radiant  and  conductive  heat  losses.  It  also  neglects 
the  heat  needed  to  bring  the  sublimating  material  from  its  initial 
temperature  to  sublimation  temperature;  both  are  fair  approximations 
for  low-temperature  sublimatcrs. 

If  we  define 


B' 


m 


p  u  C, 
e  e  h 


(h  -  h  ) 
aw  w 


(7) 


and  use  the  identity 


B_ 

B' 


T 

aw 


o 


(8) 


we  obtain  the  second  well-known  relation  of  simple  theories. 
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BOUNDARY- LAYER  MASS-TRANSFER  CHARACTERISTIC 

The  solution  of  the  boundary- layer  diffusion  equation  for  the 
transfer  of  the  sublimated  material  away  from  the  wall  is^ 


-  V 

("re  •  V 


(9) 


TMs  form  is  strictly  true  only  for  the  case  where  both  Le  and  Pr  are 
equal  to  one;  in  this  case  the  diffusion  equation  and  the  equation  for 
the  distribution  of  total  enthalpy  in  the  boundary  layer  are  identical 
and  concentration  and  energy  profiles  are  similar.  The  effect  of  Pr 
on  the  similarity  of  the  profiles  is  minor,  weaker  than  its  ef¬ 

fect  on  itself.  This  justifies  the  use  of  Eq.  (8)  without  also 
setting  Pr  *  1  in  Eq.  (1)  and  equations  derived  from  it. 

Combining  Eq.  (9)  with  Fick's  law  (see  Ref.  7),  we  obtain  an 
expression  for  the  wall  concentration  as  a  function  of  the  rate  of 
sublimation  at  the  wall: 


m  [  1  -  C  |  ■  C 
s  I  sw 


PHASE -CHANGE  KINETICS 


.  h  -  h 

H  „  aw  w 

r —  =  c  p  u  •  •  ■  -  c.  ( 10  > 

sw  h_  -  h  swHe  e  h_,  -  h  h  * 

xe  w  Te  w 


A  last  equation  couples  the  concentration  of  the  sublimated  spe¬ 
cies  in  the  boundary  layer  to  the  sublimation  phenomenon  itself.  The 
net  rate  of  exchange  of  surface  material  across  the  surface-potential 
barrier  is 


®s  "  2ttRT  <P. 


equ  s 


P  )  ■  C 


PM 


<C 


v^roTf  '  *  *9“  “ 

s  w 


C  )  (11) 


where  c  is  an  empirical  "vapo :ization  coefficient,"  P  and  C 

r  s  equ  s  equ 

are  the  equilibrium  partial  pressure  and  concentration  corresponding  to 

the  wall  temperature,  and  P  and  C  are  the  actual  partial  pressure 

8 

W  W 


and  concentration  of  the  subliming  material  immediately  over  the  surface 
in  the  boundary  layer.  For  a  two-component  mixture  one  has  the  follow¬ 
ing  auxiliary  relations  between  the  concentrations  and  the  molecular 
weights  (M  is  the  mean  molecular  weight  of  the  mixture): 


An  important  characteristic  of  the  phenomenon  is  that  there  is  a 
maximum  rate  of  escape  of  surface  atoms,  which  occurs  when  the  concen¬ 
tration  of  the  material  in  the  surrounding  gas  is  zero  and  which  de- 

(2) 

pends  only  on  the  wall  temperature.  Several  analytical  expressions 

for  the  value  of  P  that  determines  this  maximum  can  be  written 

s  equ 

down,  depending  on  the  subtlety  of  the  microscopic  model.  The  simplest 
one,  corresponding  to  the  integral  of  the  Clausius-Clapeyron  equation, 
is 


P 

8  equ 


P  exp 


(13) 


where  p  is  a  constant.  This  yields  for  the  maximum  (forward)  vapori¬ 
zation  rate  the  expression 
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The  set  of  Eqs.  (1),  (3),  (9),  and  (13),  together  with  an  auxil¬ 
iary  caloric  equation  of  state  defining  c 


h  -  h 
aw  w 

T  -  T 
aw  w 


complete  the  general  definition  of  the  problem.  We  shall  assume  in 

this  analysis  that  c  can  be  treated  as  an  a  priori  (or  iteratively) 

P  (11) 

determined  constant,  not  dependent  on  the  concentrations. 

The  following  dimensionless  parameters,  which  depend  only  on  the 

properties  of  the  surface  and/or  the  free-stream  conditions,  are  now 

defined: 


H  - 


A  H 

S  8 


M.c 
A  p 


V  -  1 


JL  £_ 

Pr 


k  *  1  - 


H1*?  -  » 


>  -  1 


5  -  -  V2T?  IT  T  o  u  C.  a  /  He 

*  tY  M  aw  eehy  x 
s  o  o  V 


The  function  §  represents  the  product  ChA/  Re  which,  at  least 

o  ^ 

for  self-similar  flows,  is  a  constant.  Otherwise  this  product  varies 
with  x  through  the  streamwise  pressure  and  the  wal l -temperature 
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gradients.  While  questions  must  be  raised  as  to  ita  behavior  near  the 
singular  point  x  ■  0,  which  are  mentioned  in  Section  IV,  at  large  Rey¬ 
nolds  numbers  it  is  certainly  a  finite,  slowly  varying  quantity. 

Eliminating  m  between  Eqs.  (2)  and  (11)  and  rearranging,  we  write 
0 


p  u  C. 
Ke  e  h 


V2nRT 

_ 3 


(21) 


By  straight  substitutions  we  transform  it  into  an  equation  in  only  one 
unknown,  which  is  B: 


For  clarity,  a  few  intermediate  steps  In  this  transformation  are 
given  in  the  following.  The  left-hand  sides  of  Eqs.  (21)  and  (22)  in¬ 
volve  simply  the  definition  of  $,  Eq.  (20).  The  first  factor  on  the 
right-hand  side  follows  from  Eqs.  (8)  and  (4),  which  yield 

B 


B'  - 


(23) 


Substituting  for  B*  from  Eq.  (7)  one  obtains 


B  - 


-If- 

'  aw 


-  1 


c  T 
p  aw 


Of 

Pr 


fr- 

\  aw 


From  Eqs.  (S)  and  (18)  one  has 


r  -  t  (i  +  tb) 

av  aw 

o 


(24) 


(25) 


which,  introduced  into  Eq.  (24)  and  solved  for  the  temperature  ratio, 
gives 


(26) 


The  first  term  in  brackets  in  Eq.  (22)  is  simply  Et.  (13)  written 
In  the  form  (see  Eqs.  (16)  to  (19)) 


P 

a  equ 


H  H 


pe  e 


(27) 


with  the  temperature  ratio  as  given  by  Eq.  (26). 

The  second  term  In  brackets  follows  from  Eqs.  (12)  and  (10). 
From  Eq.  (10)  we  derive  the  following  expression  for  C#w  In  terms 
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Thia  equation,  after  substitution  in  Eq.  (12),  lead*  to  an  expreaalon 

for  P  which  la  exactly  the  aacond  tar*  In  bracketa  of  Eq.  (22). 
iv  7 
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III.  ASYMPTOTIC  BEHAVIOR  DOWNSTREAM 


For  large  values  of  Re  Che  left-hand  side  of  Eq.  (22)  vanishes 

Sr 

and  the  equation  yields  a  non-zero  value  B  **  B^  given  by  the  solution 
of  the  factor  in  brackets  on  the  right-hand  side. 

Note  that  this  corresponds  exactly  to  stating 

P  -  P  =*  0 

s  equ  sw 

The  solution  is  identified  with  the  "equilibrium  solution"  in  which 
the  partial  pressure  of  the  vapor,  and  consequently  also  its  concen¬ 
tration  and  temperature,  are  values  corresponding  to  thermodynamic 
phase-change  equilibrium. 

In  the  present  formulation  the  solution  is  given  in  terms  of  the 

constant  p  in  the  analytical  expression  for  P_  (instead  of  speci- 

s  equ 

fying  Tw  *  T?  ^  separately  from  thermodynamic  tables).  A  convenient 
graphical  procedure  ia  obtained  defining 


(29) 


ia  terms  of  which  the  asymptotic  solution  cakes  the  form 


H 

X  ♦  1 


♦  In 


♦  k  -  PH 
X  T~k 


(30) 


This  is  plotted  in  Fig.  1  for  a  particular  value  of  f.  For  any  exter¬ 
nal  static  preasure  P#  and  parameters  des'ribing  the  properties  of 
the  sublimating  material  (p  and  H),  one  obtains  a  unique  value  for  the 
well -temperature  function  Xm.  With  this  value  of  Xm  and  further 


*The  singularity  B  •  0  is  uninteresting.  It  corresponds  simply 
to  flow  without  sublimation. 


k  =  0{Pf - 1) 

k  *  -0 1797  (Pr  =  0.7,  M>8) 
r*  0.287  (Qir) 


Asymptote  X-,  =  -k  or  T  *  T 
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parameters  pertaining  to  the  reference  flow  and  the  sublimating  mate¬ 
rial  (a/Pri-),  Fig.  2,  which  is  a  plot  of  Eq,  (29),  yields  the  value  of 
the  bloving  parameter  B0. 

Figure  3  is  an  auxiliary  figure  giving  the  variation  of  Q'/Prr. 
This  can  be  written  as 


or 

Prt 


(31' 


The  two  proportionality  constants  or  and  g  (see  Eqs.  (4)  and  (5))  ap¬ 
pear  as  a  ratio.  Therefore,  their  dependence  on  the  principal  varia¬ 
bles  of  the  binary  bound ary -layer  problem,  which  are  M  /M  ,  Pr,  and 

8  A 

pressure  and  wall -temperature  gradient,  is  minimized.  If  one  makes 
the  hypothesis  that  both  vary  in  the  same  fashion  with  these  parame¬ 
ters,  which  seems  probable,  then  the  combination  a/PrT  Is  to  a  very 
good  approximation  a  function  only  of  the  free  stream. 


40 


Fig.  2— Asymptotic  blowing  parametei 


9 


8- 


Fig. 3 -Variation  of  some  Mach -number -dependent  parameters (tt-0) 
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IV.  BEHAVIOR  NEAR  THE  STAGNATION  POINT 
Equation  (31)  has  the  form 

-  (32) 

The  parameter  §  appearing  on  the  left-hand  side  of  Eq.  (32)  rep¬ 
resents  the  product  C^^Re~ .  The  assumption  is  made  that  this  product 

"o  x 

is  a  constant,  as  it  is  for  self-similar  boundary-layer  solutions. 

Now,  since  m  ■  Xq  is  bounded  according  to  the  surface-evaporation 
s 

equation,  the  assumption  implies  chat  the  wall  temperature  tends  to 
the  recovery  temperature  towards  the  leading  edge  (Re^  -0).  It  fol¬ 
lows  that  the  wall  temperature  is  variable,  and  C  *  constant  can 

no  * 

only  be  good  in  the  sense  of  "local  similarity." 

It  is  worth  noting  that  the  above  implies  two  statements  of  un¬ 
equal  importance  to  the  present  analysis.  The  more  Important  one  con¬ 
cerns  the  behavior  of  C,  ;  that  is,  that  C.  grows  without  bounds  to- 

o  "o 

wards  the  leading  edge.  This  leads  to  the  result  that  B  ~  mg/Ch  tends 
to  aero  there  and,  thus,  that  B  varies  between  zero  and  B  over  the  en- 

OP 

tire  region  of  interest.  The  second  and  less  important  one  concerns 

the  numerical  accuracy  of  the  assumption  C.  "\f  Re  *  constant.  It  is 

"o  x 

undoubtedly  poor  very  near  the  origin  where  the  wal 1 -temperature  gradi¬ 
ents  are  large,  but  it  is  probably  satisfactory  in  the  region  of  as¬ 
ymptotic  approach  to  the  near-equilibrium  solution  downstream  of  the 
leading  edge. 

It  is  not  possible  to  discuss  conclusively  the  difficult  problem 
of  the  singularity  at  the  leading  edge.  It  must  be  accepted  on  the 
basis  of  heuristic  arguments  and  the  analogy  with  the  behavior  of  or¬ 
dinary  boundary  layers  at  the  leading  edge,  which  involves  similar 
difficulties.  Physically,  the  behavior  outlined  in  what  preceded  is 
quite  reasonable.  Moving  upstream  towards  the  leading  edge,  the  heat 
flux  to  the  wall  increases,  and  the  wall  temperature  must  rise  to 
permit  an  increased  rate  of  sublimation. 
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Returning  to  Eq.  (32),  the  right-hand  side  is  expressible  in  the 
form  of  a  series  in  the  interval  0  <  B  <  B^.  If  the  asymptotic  blow¬ 
ing  parameter  is  reasonably  small,  which  is  already  implied  by  the 
use  of  the  linear-blockage  equation,  Eq.  (4),  a  limited  number  of  terms 
of  the  series  will  represent  the  behavior  of  the  function  throughout 
the  range. 

The  series  is 


urn.  m  i  +  noi+i 

B  f  (0)  2  f  (0) 


(33, 


with  the  following  expressions  for  the  function  and  its  derivative 
evaluated  at  the  origin: 


f(0) 


(34) 


iL m 

f(0) 


(35) 


The  parameter  a  has  the  form 


(HO2  pe~H  +  k(k  -  HO  fHpe~H 
pe*H  ♦  k(k  -  2H0  pe'H  ♦  kP  (PH  -  k) 


(36) 


and  it  was  defined  so  as  to  become  one  when  k  *  0  (Prandtl  number  uni¬ 
ty). 

It  is  interesting  to  record  the  initial  behavior  of  the  pertinent 
physical  parameters  of  the  problem.  First-order  expansions  arc  given 
below: 
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sw 


k 


k 


r~H 


M 

s 


1 


P  k 

P  «  — £ - 

sw  k  -  rn 


p 

S  Cqit 


(37) 


(38) 


(39) 


(40) 


(41) 
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V.  SCALE  OF  THE  LF*T)ING-EDGE  EFFECT 

The  first  two  terms  of  the  series  expansion,  Eq.  (33),  vary  as 
Re^  -*  ®  to  the  asymptotic  limit  of  B  ■  f(0)/f'(0).  This  limit  is  not 
correct  because  it  does  not  equal  B^  unless  tends  to  zero.  However, 
we  know  independently  from  the  full  solution  at  the  asymptotic  limit, 
Eqs.  (29)  and  (30),  or  Figs.  I  and  2. 

The  argument  suggests  that  if  one  replaces  in  the  two-term  ex¬ 
pansion  Eq.  (33) 


Lm  .-L 

f(0) 


(42) 


one  obtains  a  good  approximation  to  the  behavior  of  B  fir  all  Re^  and 
arbitrary  values  of  B  ,  provided  is  sufficiently  small  to  justify 
the  linear  expression  for  the  blockage  factor,  Eq.  (4). 

Accordingly,  an  approximate  equation  for  B  is  as  follows: 

(43a) 


B_ 

P. 

•r- 


(43b) 


where  u  is  i  dimensionless  scale  factor 


(44) 

The  corresponding  approximation  for  the  variation  of  mass-addition  rate 
is,  from  Eq.  (45) 
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m 

s 


-R 

X 

V  Re  +  J~7 

x 


(4>) 


For  pointed  bodies  the  trends  under  study  are  demonstrated  more 
emphatically  if  one  defines  a  length  by 


*  u 

e  e 


(46) 


which  is  explicitly  given  as  a  function  of  reference -f low  and  subi 

'k 

mating-n ater i al  parameters  by 


Equation;  (43)  and  (45)  can  tin  '  be  written 


R 


♦  1 


(48) 


ILi £1 


(49) 


This  i  ;u  t  ion  shows  iww,  for  flow  i  roo  a  pointed  leading  edge.  *a 
starts  i  ra-.  its  -i.!*  imu-  v  iu*  dictated  by  kinetic  considerations  and 
i!vcr>  as*  i  to  Its  asymptotic  variation  downstream,  which  is  proportion.! 


hot «  that  C  \  Re 
h  x 

o 


is  a  constant  with  x. 
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to  the  inverse  square  root  of  the  streaatwise  distance. 

Similar  explicit  expressions  for  the  variation  of  the  other 
parameters  are  obtained  by  substituting  B  from  Eq.  (48)  into  the 
first-term  expansions,  Eqs.  (37)  through  (40).  We  see  that  the  con¬ 
centration  and  the  partial  pressure  of  the  sublimating  species  reach 
constant  values  eit  the  origin,  the  magnitude  of  which  depends  only 
on  k  (the  Prandtl  number)  and  which  are  zero  for  the  particular  case 
of  Pr  =  1,  The  wall  temperature  tends  to  the  adiabatic  recovery  tem¬ 
perature  for  zero  blowing  at  the  origin;  that  is,  to  the  stagnation 
temperature  if  Pr  =  1.  The  blowing  parameter  itself  tends  to  zero  at 
the  origin  of  the  boundary  layer. 

The  first-term  expression,  Eq.  (48),  is  plotted  in  comparison 
with  a  numerical  solution  of  the  full  Eq .  (22)  in  Fig.  4  for  a  partic¬ 
ular  set  of  physical  parameters. 

It  is  of  interest  to  consider  further  the  parameter  A,  which  de¬ 
termines  the  scale  of  the  adjustment  to  the  asymptotic  solution.  Fig¬ 
ure  5  shews  the  variation  of  this  quantity  in  the  form  (see  Eqs.  (44) 
and  34)) 

Vap~  B 

T*  *  T - 2 - £  (50) 

p" exp  ("H)  “  F^Th 

e 

The  principal  independent  parameters  in  Eq.  (50)  are  p/Fg  and  H 
(note  k  *  0  for  Pr  =  1).  These  also  determine  B^  uniquely  (see  Pigs. 

1  and  2),  provided  that  the  secondary  parameters  or/PrT  and  t  can  be 
considered  constant,  which  is  indeed  very  closely  true  above  a  Mach 
number  of  8  (see  Fig.  3). 

Equation  (50)  is  plotted  in  Fig.  5,  and  values  of:  B  are  also 

shown  to  bring  out  the  physical  problem.  The  domain  of  practical  in- 

teresc  is  at  high  values  of  p/'P  and  relatively  large  values  of  B  , 

e  » 

The  maximum  of  these  curves,  which  occurs  when  B  tends  to  zero,  is 

1  iX)  3 

not  interesting  in  itself,  but  it  serves  to  demonstrate  the  influence 
of  another  parameter  of  Importance  which  is  the  mass  ratio  M  /M  .  In- 

S  A 

deed,  in  this  region  a  simplified  explicit  solution  can  he  written. 


Fig.  4  — Comparison  of  full  solution  with  the  approximate 
modified  first- term  expansion 
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Fig.  5 -Typical  variation  of  the  leaoing  -  euge  scale  parameter, 
showing  lines  of  -  constant 
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Low  values  of  p /P  correspond  to  small  values  of  (see  Eq.  (30)). 
Since  k  is  normally  near  zero,  the  approximate  form  of  Eq.  (30)  is 


-(x  +  k)  -  IDl .-.i0p,e 

x  00  '  P 

e 


(51) 


Moreover,  when  X  tends  to  zero,  B  tends  to  zero  also.  One  can 

*  <B  ’  CD 

neglect  the  quadratic  term  in  Eq.  (29)  which  becomes 


a 

Pr 


B 

00 


cr 

PrT 


X 


00 


Di  _A 

T  M 


s 


(52) 


Solving  for  B  in  terms  of  p/P  from  Eqs .  (51)  and  (52),  the  following 

°  oo  C 

approximate  form  of  Eq.  (50)  is  obtained: 


(H/M  JH) 

— V5 (53) 

(VV™ 


'J*. 


M. 


'max 


1  - 


This  relation  is  plotted  in  Fig.  6. 

Finally,  the  results  can  be  put  into  another  form  which  has  a 
direct  physical  significance.  Recasting  the  equations  in  terms  of 
the  length  scale  6  of  the  leading-edge  effect  by  combining  Eqs.  (47) 
and  (50)  and  using  the  equation  of  state,  one  obtains 
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2  S 
PC  ST 

L  MAj 

M  _L. 
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v  e 

0 

(54) 

2 

V  ho  ') 

r*T  t 

awo 

T  Me 
L  e  J 

n  a 
e  e 

■>  m 

P 

e 

§ 

«■ 

p  o 
—  B2 

p  ® 


exp  (2H) 


(54) 
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Fig.  6 -Behavior  of  the  maximum  of  the  scale  length 
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This  is  plotted  in  Fig.  7.  Note  that  the  three  bracketed  terms  of 
the  denominator  on  the  left-hand  side  depend  only  on  the  basic  flow 
field,  the  Mach  number,  and  the  free-stream  static  temperature,  re¬ 
spectively.  The  bracketed  term  in  the  numerator  depends  only  on  the 
properties  of  the  sublimating  material.  The  entire  factor  on  the 
left-hand  side  is  independent  of  pressure.  Consequently,  Fig.  7  shows 
that  during  the  initial  re-entry  from  space  (p/Pg  *  00 )  during  which 
the  Mach  number  is  roughly  constant  (therefore  H  ■  constant  and 
C  *  constant)  ,  6  increases  to  a  maximum  at  some  altitude  fixed  by  the 
value  of  the  parameter  p  and  then  decreases. 


Fig.  7 -Typical  variation  of  the  rate -controlled  flov 


length,  showing  lines  of  -constant 


-31- 


VI.  DISCUSSION 


In  order  to  show  the  magnitudes  involved,  consider  the  following 

(9) 

estimates:  Let  the  model  be  a  slender  cone  for  which 


ch -/rT  -  0.512 

o 


Assume  ^  and  a^  constant  at  their  value  in  the  tropopause  (^a^  - 
4  x  10*4  lb/ft).  It  follows  that 


aw 
_ _ o 

T  Me 


-  a 
e  e 


1.3 M*  x  10’4  lb/ ft 


(55) 


The  constant  p  can  be  obtained  from  vapor-pressure  data.  Reference  12 

(pp.  1751  -  1755)  gives  a  table  for  selected  organic  and  inorganic 

substances  from  which  it  is  evident  that  a  representative  value  is 
9  2  * 

p  *  10  lb/ft  .  (The  value  p  varies  around  this  average  by  one  order 
of  magnitude  for  almost  all  the  substances  listed.)  The  vaporization 
coefficient  C  is  poorly  known.  For  solid  sublimators  it  is  doubtful 
that  it  should  exceed  0.1,  '  and  it  can  be  less  than  that  by  two  or¬ 

ders  of  magnitude.  Finally,  the  ratio  M  / M.  can  be  taken  as  unity  for 

S  A 

the  purpose  of  this  estimate. 

It  follows  that,  quite  generally, 


10 


-11 


M> 


10  7  ft 


(56) 


The  magnitude  of  is  seen  to  depend  very  strongly  on  the  value 
of  K.  Since  it  increases  rapidly  with  H,  let  us  illustrate  the  prob¬ 
lem  for  a  large  value  of  H:  Data  typical  of  graphite  (*  •  25,000  Btu/lb, 

M  -  12)  at  a  flight  Mach  number  of  15  yield  approximately 
s 

*  b 

In  the  notation  of  the  reference,  p  ■  10  where  b  is  tabulated. 
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16.7 


9  2 

The  value  of  p  »  10  lb/ ft  used  in  the  preceding  estimates  represents 
graphite  quite  well--well  enough  considering  the  uncertainty  in  the 
evaporation  coefficient  «.  Since  6  increases  with  pressure,  let  us 
consider  relatively  low  altitudes  consistent  with  the  assumed  Mach 
number  and  high-speed  re-entries,  say  30,000  ft.  Consider  a  20-deg 
cone  for  which  the  surface  pressure  is  approximately  300  P^ .  It  fol¬ 
lows  that 


10*  3  <  5  <  10’1  in. 

Had  an  altitude  of  25,000  ft  been  is?  rr-od ,  t  would  be  larger  by  an  or¬ 
der  of  magnitude.  Had  data  typical  of  teflon  been  used  (*.  ~~  1000  Btu/lb, 

M  ■  100),  H  would  have  been  in  the  order  of  5,  and  5  would  be  entirely 
s 

negligible. 

It  is  interesting  to  note  that  the  maximum  .  f>r  a  given  H  occurs 
at  values  of  p/P^  which  correspond  to  considerably  higher  than  1 
atm.  While  pressures  ab^ve  atmospheric  can  occur,  because  p^  and 
are  values  outside  the  boundary  layer  behind  the  leading  shock,  they 
would  not  normally  be  as  igh  as  indicated  for  i  max. 

The  rate-controlled  region  can  be  taken  to  extend  over  a  distance 
from  the  origin  equal  to  100  5;  that  is,  for  points  on  the  surface 
lying  beyond  this  limit  the  error  incurred  by  using  the  equilibrium 
solution  for  B,  T^,  etc.,  is  less  than  10  per  cent  (see,  for  ex¬ 

ample,  Fig.  4  or  Eq.  (48)).  The  sketch  on  the  following  page  shows 
the  behavior  of  the  rate  of  sublimation  indicated  by  both  the  equi¬ 
librium  approximation  snd  the  full  solution  and  illustrates  the 
present  argument. 

The  preceding  estimates  of  the  scale  of  the  leading-edge  effect 
indicate  clearly  chat,  in  comparison  wetb  reasonable  man-saadc  re-entry- 
vehicle  sires,  the  region  of  transitional  sublimation  is  very  snail, 
indeed  normally  negligible.  However,  f  is  ar.  independent  parameter 


Equilibrium  solution,  mg  — 

V* 


Pr*s«nt  solut'on,  mt  — 


7F 


Equilibrium  approximation 
valid  withm  10  per  cent 


Rote-controlled 
region*"*-- ! 


x  =  1008 


x  (distance  from  the  ieoding  edge) 


which  does  not  scale  with  the  body  geometry.  Th  ju-imation  of  smal 


models  m  wind  tunnels  or  ireteorites  in  the  "t^sw.-cre  aay  be  total  lv 


dominated  bv  transitional  effects. 


This  scaling  problem  can  best  be  demonstrated  by  using  a  concrete 
example.  Consider  the  w.nd- tunnel  experiments  reported  in  Ref.  3  using 
camphor  in  a  Mach  3  wind  tunnel.  For  camphor  the  materials  data  at 
actual  lest  conditions  can  be  given  quite  accurately,  except  for  C: 


9  '> 

330  j/gm,  -  132,  p  ■  1.74  x  10  Ib'ft*.  Estimate  «  .»*•  0  1  as 


be  for* 


With  the  tunnel  recov  rv  temperature  at  3S0°K  and  the  static 
pressure  of  10  m  Hg,  one  finds  M  -  17.2  »j.d  p/p  •  fc  ;  k  10#: 


2 . h  x  10  in 


Considering  that  the  transitional  sublimation  region  extends  to  a  dir  - 
tance  of  the  order  10O.  (see  Fig.  4),  one  concludes  that  the  ent  i  r<= 
v l nd - 1 ur.ne  1  mod e  1  is  affected  by  transitional  effects  u;  ier  these  test 


c  ond  v  t  i  or  s 
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One  other  aspect  of  the  aerodynami  i  of  subliming  bodies  can  be 
affected  by  the  present  resales  in  an  important  manner.  It  is  the 
problem  of  surface  recession  and  terminal  shape.  Briefly,  the  shape 
of  a  slender  body  ^  at  any  instant  t  is  given  by  the  solution 
(with  proper  boundary  conditions)  of 


at 


m 


/■J 


.< 


* 


r~ 


x  -  V 


(t) 


-r  V 


where  <  is  o  constant,  X*  the  distance  to  the  point  (x,y)  from  the 

leading  edge,  which  is  itself  receding  relative  to  fixed  coordinates 

at  an  unknown  rate  Z .  v . 

it; 

The  nature  of  this  problem  is  such  that 


Lim  y  f 


6-0 


^  > 


*  y(x,t,6=0) 


and  the  analytical  solution  for  the  terminal  shape  of  pointed  bodies 

using  the  riear-c  ;ui'  '  rium  form  of  the  local  sublimation  rate  is  not 
* 

correct . 


* 

A  study  of  this  problem  is  in  progress  at  The  RAND  Corporation. 
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